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Abstract 

The high temperature phase of the three dimensional random field Ising model is studied 
using replica symmetry breaking framework. It is found that, above the ferromagnetic transi- 
tion temperature Tf , there appears a glassy phase at intermediate temperatures Tj < T < Tb 
while the usual paramagnetic phase exists for T > Tb only. Correlation length at Tb is com- 
puted and found to be compatible with previous numerical results. 



Although a great deal of work has been devoted to the understanding of the random field 
Ising model (RFIM) |l|], some aspects still need to be cleared up. It is now well-known that, in 
dimension D = 3, long range order is present at sufficiently low temperature and weak random 
fields with non trivial critical exponents ||, the upper critical dimension of the RFIM being 
D = 6. Nevertheless, perturbation theory leads to dimensional reduction (critical exponents 
are incorrectly predicted to be equal to those of the corresponding pure model in dimension 
D — 2) H and therefore does not succeed in describing the critical behaviour of the RFIM. The 
reason of this failure presumably stems from the very complicated energy landscape due to the 
quenched disorder, and more precisely, from the existence of a huge number of local minima 
of the free energy in the space of local magnetisations that usual perturbative expansions do 
not take into account Q. Numerical simulations and resolutions of the mean- field equations 
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corroborate this picture JH, |f§. Above the ferromagnetic transition temperature TV, there seems 
to appear an intermediate "glassy" regime for TV < T < Tb where many solutions of the local 
magnetisations mean-field equations coexist, while only one of them subsists in the paramagnetic 
phase T > Tb- From the theoretical point of view, it was suggested that the techniques of 
replica symmetry breaking (RSB), which proved to be successful in the mean- field theory of 
spin glasses J7|] where such complicated free energy landscapes arise, could also be applied to 
the RFIM ||. Experiments made on diluted anti-ferromagnets also found an irreversibility line 
above the critical temperature where the anti-ferromagnetic order appears Q. Recently, Mezard 
and Young, referred to in the following as M— Y, proposed a variational approach of the RFIM 
|l0| based on a self-consistent expansion in 1/N (where N is the number of spin components) 
due to Bray [O]. They found that replica symmetry, which gives back dimensional reduction, 
must be broken at the ferromagnetic transition T = TV and that RSB solution leads to sensible 
results for the critical exponents in agreement with already known results |i"of . 

In this paper, using M— Y's framework, we concentrate upon the non ferromagnetic regime 
(i.e. T > TV). We find that there exist indeed two different phases : a paramagnetic phase 
at high temperatures T > and a glassy phase at intermediate temperatures Tf < T < T^. 
The value of the correlation length at T = where the RSB transition occurs is computed and 
compared to predictions obtained from numerical resolution of mean-field equations ||. 

The model we consider is a ^-component version of the RFIM on a three-dimensional lattice 
including L 3 spins = ($L where i = (ii, 12, is) and < ii, 13 < L — 1, 

h) = \ £ (*< - 2 + I £ (* -) 2 + ^ E ((*i) 2 ) 2 " E (1) 

i i i 

where is a quenched random field, the distribution of which is Gaussian, uncorrelated at 
differents sites, with mean = and variance h^.h x - = A <5^ A 8a a. Following the standard 
procedure (7J, we introduce n replicas of the spins <& a , a = l...n, and average over the quenched 
disorder h to obtain the effective Hamiltonian 

n{* a }) = \ _E (*?-*5) 2 + iE(*l) 2 -f E *t-*i + ^E((>i) 2 ) 2 (2) 
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The replica correlation functions (^i^^'k^ = ^ ab (^) w ith the effective Hamiltonian (|2|) are 
related to the disconnected and connected correlation functions, (^fc)h-(& X fc)h = ^ X Gdi s {k) 
and {$k-^ X _k)h - (®k^ h '(®-k) h = ^ G con {k), where (.)h denotes the average over the Gibbs 
measure induced by (jl|) |i0| . For the replica symmetric assumption G ab (k) = G(k) 5 ab + G(k), 
the correspondance is simply Gdi s {k) = G(k) and G con (k) = G(k). 

Using Bray's self-consistent screening approximation [fLl| l which is exact to order 1/N, one 
finds that the propagators G ab (k) are given by the saddle-point of the free energy 

F({G«\k)}) = E^-2Ecos(^)+r + ^3E Gaa (9)-AjG-(fc) 

/ \ aa 1 / \ aa 

- 2 A E G ab (k) - E log G(k) + - E log [1 + n(fe)] (3) 
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where 



n a \k) = ±Y, Gab ( k -<i) Gab (<i) 



(4) 



Within the replica symmetric hypothesis, the correlation functions are therefore solutions of the 
following set of 2L 3 implicit equations 
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In order to determine where the transition to RSB occurs, we have studied the local stability 
of the free energy T(G ab ) around the symmetric saddle-point (G,G). Repeating de Almeida— 
Thouless's calculations |l2| and taking into account the fc-dependence of the order parameters 
G and G, we have found that the replica symmetric solution is locally stable if and only if the 
lowest eigenvalue A of the matrix 
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is strictly positive. 

When N — > oo, Bray's partial resummation reduces to the Hartree-Fock approximation. 
G(k) is thus equal to the bare propagator with a renormalised squared mass m 2 solution of 
the gap equation (||, |6|). From the expression of the AT matrix (P), one obtains A = m 4 > 0. 
The replica symmetric solution is therefore always stable. As soon as iV becomes finite, the 
corrections appearing in (||) may lead to instabilities. In our three-dimensional system, however, 
the self-consistent screening approximation induces no ferromagnetic transition for large iV (m 
never vanishes for finite bare temperatures r) |10| 11]. Hereafter, we choose N = 1 (Ising case) 
which allows for the existence of long range order at finite r < 0. 

For every size L of the lattice, we fix a value of r and solve for the propagators G, G by an 
iteration of equations (||, ||). Using rotational and translational symmetries, only G(k) and G(fc) 
with < ki < k2 < k% < Int(^) are to be found. Once a fixed point is reached, we estimate the 



mass m and the correlation length £ 



from the low-momentum behaviour of G 



G(0,0,0)^ — 



(10) 



Expression (|l(]) is exact for N = oo ^ 2 '. Moreover, from perturbation theory which is thought to 
be correct above the AT transition, we expect G(k) to have a single pole. A similar calculation 
gives m and £ = assuming that G(k) has a double pole ( 3 ). The lowest eigenvalue A is then 
computed by diagonalising the AT matrix. This highly time consuming task may be simplified by 
observing that A4 is invariant under the three symmetry operators Sd : — > L — kd- In the base 
of their eigenvectors, A4 reduces to eight diagonal blocs of size roughly equal to^-r^ x (l 1 ) • 
On can check that the eigenvector corresponding to A belongs to the "physical" subspace, i.e. 
the one spanned by eigenvectors of Sd of eigenvalues +1. The process is repeated until the value 
Ti, of the bare temperature r where A vanishes and the corresponding correlation lengths £l and 
£l are bracketed with a sufficient precision. The final uncertainties on £l and £l are lower than 
±5.10 -4 for lattice sizes running from L = 2 up to L = 20. 

The numerical values of the correlations lengths at the AT transition are displayed Fig. 1. 
Although it seems difficult to extrapolate to L — * oo, reliable information on the thermodynam- 
ical limit may be obtained since the correlations lengths are relatively small as compared to the 
lattice size (£l < £x < L/3 for L = 20). From finite size effect theory [l3|], we expect indeed 
that, if the mass tHl converges to a finite value > at the thermodynamical limit, then its 
asymptotic behaviour obeys 



rhL-rhoo ~ C . e- L/ ^ L + O (e^ 2L/ ^) (11) 



where C is a constant (the same identity holds for m and £l with a different constant C). Fig. 2 
shows the dependence of rhi and mj, upon e~ L ^ L and e~ L ^ L respectively. The linear law ([Tl| ) 
is very well verified with proportionality factors of order one (C ~ 1.83, C ~ 1.37) W. Linear 
extrapolations to L — > oo provide the values of the correlation lengths at the thermodynamical 
limit 

loo =i Coo =i 7.7 ±0.2 (12) 

The equality between the correlation lengths defined from the disconnected and the connected 
correlation functions is a self-consistent check of our analytical and numerical results. It is 
indeed predicted by perturbation theory [^, ||] and therefore holds for high temperatures down 
to the RSB transition. 

In this letter, we have argued that the non- ferromagnetic phase of the three-dimensional 
random field Ising model is composed of a paramagnetic phase at high temperatures and a 
spin-glass phase at lower temperatures. The onset of this glassy phase therefore occurs at a 
finite correlation length for both correlation and susceptibility functions which was found to be 



2 For e.g. N = 4, long range order is absent and both mass m and eigenvalue A are always non zero. However, 
from numerical resolution of the saddle-point equations (|^, we found that using definition (|l(j) of fh, the 
relation A ~ m 4 is roughly correct. 

3 This stems from perturbation theory which is thought to be valid above Tb ^) and can be checked up when 
TV — *• oo on ^ ^). For the Ising system, the one dimensional correlations functions in the real space, g(x) and 
g(x) (which are respectively the Fourier transforms of G(fc,0,0) and G(k, 0, 0)) can be fitted with a very good 
agreement by g(x) ~ acosh(x/f) + b and g(x) ~ a(l + x/£) cosh(a;/£) + b ||. The discrepancy between these 
values of the correlation lengths and the ones defined in the text seems to vanish for increasing lattice sizes. 

4 Although scaling relation (|l"l| ) should also be valid for the correlation length itself, i.e. foo — £x ~ D . e~ L ^ L , 
one sees from D — —C. ^ ~ —100 that it cannot be directly observed with lattice sizes lower than L — 20. 



in the range 7.5 < £ < 8. Although such a result might be an artefact due to the l/N approach 
used here, it is in qualitative agreement with previous numerical studies which found that the 
mean field equations of the RFIM begin to have more than one solution, and thus that the 
perturbative approach ceases to be correct, for £ > 4.5. Such a behaviour is to be expected in 
the low temperature phase too. Since at very low T, only the two states where all spins are 
aligned along the same direction remain, replica symmetry has to be restored at a temperature 
T s with < T s < Tj. It would be interesting to extend the calculation we have presented here 
to verify explicitely this conjecture. 
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Figure Captions 

Figure 1 : the correlation lengths £l (empty dots) and £l (full dots) at the onset on the spin 
glass phase for different lattice sizes L = 2 up to L = 20 (the total number of spins is L 3 ). 

Figure 2 : the masses fhi = j- (empty dots) and mi = |- (full dots) plotted vs. the "finite 
size" factors exp(— L/£l) and exp(— L/^l) respectively. The dashed lines are the best linear fits 
from the last twelve points (sizes L = 9 to L = 20). 




Figure 1: correlation lengths £l (empty dots) and £l (full dots) at the RSB transition for 
differents lattice sizes L. 
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Figure 2: masses tHl = j- (empty dots) and itil = j- (full dots) vs. e~ L ^ L and e~ L ^ L 
respectively. The dashed lines are the best linear fits from the last twelve points (sizes L = 9 to 
L = 20). 



